function [U_11, U_12, U_21, U_22,...
          V_11, V_12, V_21, V_22,...
          Phi_11, Phi_12, Phi_21, Phi_22,...
          D, AvU21, AvV21]=load_function(fp0)

filepath = fullfile( fp0, 'U_intervals.mat' );
A=load(filepath);
U_11=A.U_11;
U_12=A.U_12;
U_21=A.U_21;
U_22=A.U_22;

filepath = fullfile( fp0, 'V_intervals.mat' );
A=load(filepath);
V_11=A.V_11;
V_12=A.V_12;
V_21=A.V_21;
V_22=A.V_22;


filepath = fullfile( fp0, 'Phi_intervals.mat' );
A=load(filepath);
Phi_11=A.Phi_11;
Phi_12=A.Phi_12;
Phi_21=A.Phi_21;
Phi_22=A.Phi_22;

filepath = fullfile( fp0, 'D_interval.mat' );
A=load(filepath);
D=A.D_interval;

filepath = fullfile( fp0, 'Av_interval.mat' );
A=load(filepath);
AvU21=A.AvU21_interval;
AvV21=A.AvV21_interval;
end